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Abstract 

The distribution of price returns for a class of uncorrelated diffusive dynamics is considered. The 
basic assumptions are (1) that there is a "consensus" value associated with a stock, and (2) that 
the rate of diffusion depends on the deviation of the stock price from the consensus value. We 
find an analytical expression for the distribution of returns in terms of the diffusion rate, when 
the consensus value is assumed to be fixed in time. The analytical solution is shown to match 
computed histograms in two simple cases. Differences that result when the consensus value is 
allowed to change with time are presented qualitative explanations. 
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I. INTRODUCTION 



Statistical analysis of a wide range of physical and other processes have shown the presence 
of non-Gaussian distributions. The_Y, iiiclude temperature fluctuations in hard turbulence Q| , 



diffusion in inhomogeneous media 



One 



, 0| , and price variations in financial markets 
common characteristic in these distributions is the presence of exponential or power-law tails, 
signifying a more frequent occurrence of large deviations than expected from a collection 
of independent, identically distributed events |6j]. The width of the distributions have also 
been shown to scale as the time interval during which the fluctuations occur. Based on these 
properties, it has been proposed that Levy distributions be used to describe fluctuations in 
the underlying processes p,!^. In this paper, we propose an alternative explanation for the 
non-Gaussian distributions, namely a non-uniform diffusion rate 0,^^. The discussion here 
is based on fluctuations in financial markets. 

Consider a stock whose price at time t is given by S{t). Most financial market analyses are 
conducted in terms of the "return" of stock, x{t) = ln[S'(t)/5'o], where 5*0 is a reference price. 
Studies have shown that successive fluctuations in return are uncorrected 14, 111 I . Most 

n 

theoretical analysis assume further that successive events are statistically independent p. 
The traditional theoretical model uses a simple discrete random walk in the logarithm of 
price to describe the dynamics of an asset in a market. This model results in a normal, or 
Gaussian distribution of price returns |12]. On the other hand, empirical data clearly show 
that the price return distribution of real markets deviates significant Iv from a Gaussian, 



especially far from the mean 



El 



particular, some detailed studies j4 , 



(For recent reviews see 



14 
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18|.) In 



have found that the tails of the distribution have 



3,0 



an asymptotic power law decay, typically with exponents larger than 2, while others 
have found that the tails of the distribution are better described by an exponential decay. 

It has been conjectured that the non-normality observed in real financial markets can 
be explained by assuming that although successive events may be uncorrelated, the rate 
of trading depends on the price of the stock In particular, it is assumed that there 

is a "consensus" value x{t) of the return, at which point the stochastic fluctuations are 
minimum, and the level of the fluctuations increase as the return moves away from x{t). In 
other words, the rate of transactions depend on the return. In particular, a large change in 
price is likely to be followed by large fluctuations. 



2 



Here, we explore the implications of the conjecture by studying variations of the tradi- 
tional model. Like the traditional model, we use an uncorrelated random walk to represent 
the dynamics of price returns. However, unlike the former, where the rate of trading is inde- 
pendent of price, our model assumes a rate which increases the deviation of the stock price 
from S{t). It will be demonstrated that this modification in the dynamics can reproduce 
the range of non-Gaussian behavior of the price return distributions observed empirically in 
real markets. 

If an uncorrelated random walk models fluctuations of price returns, then the continuous 
dynamics is a diffusive process. In models considered here, the diffusion coefficient D{x, t) is 
assumed to depend on the time elapsed t and the return x{t). It is assumed further that the 
dependence on x(t) is a function of its deviation from the consensus return x{t) = \ii[S{t) / Sq]. 

The remainder of the paper is structured as follows. First, exact analytical results are 
presented for the distribution of price returns, valid for an arbitrary functional form of the 
diffusion coefficient, when S is fixed in time. The analytical solutions for the distributions 
are obtained by solving the Fokker-Planck equation for the corresponding diffusive process. 
Then, two specific functional forms for D{x,t) are considered. The first form is a piecewise 
linear function oi u = x/^/t. While the second is a quadratic function of u. The form of 
the distributions are obtained analytically. Results from simulations of the corresponding 
discrete random walks are also presented and shown to be consistent with the analytical 
forms. Subsequently, we provide results from computations for a model with time dependent 
consensus value. In conclusion, we summarize the results and discuss their importance for 
understanding the behavior of real markets. 



II. EXACT ANALYTICAL SOLUTION FOR THE RETURN DISTRIBUTION 



To obtain an analytical expression for the price return distribution W{x, t) of the diffusive 



processes we are considering, note that it satisfies the Fokker-Planck equation 

dW ^, ,dW 1 52 ^^^^^^ 
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(1) 



dt dx 2 dx"^ 

where D = D{x,t) is the diffusion coefficient 2^, and R{t) is a (time-dependent) drift 
rate For simplicity, we assume R{t) = for the rest of this analysis. However, the 

case of non-zero R{t) can also be treated using a simple coordinate transformation x' = 



3 



x-J^Rit') dt'. 

A normalizable solution to Eq. consistent with empirical investigations of financial 
markets [J, 0, |3| , can be found by assuming that the distribution of returns has the scaling 
form 



Wix,t) = -Fiu). 



(2) 



Here u = x/f, and rj is the self-similarity exponent P|. We also assume that the diffusion 
rate is a function of u. This scaling hypothesis leads a unique value for rj, which can be seen 
by noting that, using it, Eq. ((H) becomes 



Fiu) 



''tuF\u) = \^^{DF)"{u). 



(3) 



Consequently = 1/2, a value which consistent with conclusions from empirical studies of 
real markets. Then Eq. simplifies to 



[D(n)F(n)]"+[MF(n)]' = 0, 



(4) 



which can be integrated to 



[D{u)F{u)^ + uF{u) = Const. 



(5) 



If D{u) is symmetric about u = and the diffusion process starts at the origin, then F{u) 
will also be symmetric about u = 0. Under these conditions, both terms in the LHS of Eq. 
(0) are anti-symmetric about u = 0, and thus Const = 0. Therefore, 



D{u)F{uy = -[u + D{uy]F{u) 



(6) 



which has a general solution of the form 



exp 



u 



Diu) 



du 



(7) 



As an example, consider a constant diffusion coefficient D{x,t) = Dq. In this case, Eq. 
([7j) gives a solution of the form 

1 



F{u) = Co exp 



2Dn 



-u 



(8) 



which is the well-known result for the traditional model of distribution of returns. 
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III. STATIC CONSENSUS PRICE S 



In this section we present two models that use a non-constant diffusion coefficient D{u) 
and take the consensus value S{t) to be fixed, without loss of generahty, at S{t) = Sq. For 
both cases we present numerical and analytical results. 

A. Piecewise Diffusion 

The first form of D{u) we consider is a piecewise linear function 

Diu)=Doil + e\u\), (9) 

where Dq and e are constant parameters. It should be noted that Dq can be eliminated by 
a suitable rescaling of time. The exact solution to Eq. dT}, obtained using Eq. ((7j), is 



F{u) = Cq exp 



\u\ 



e\u\ + ir-\ (10) 



Doe_ 

where a = l/(Doe^), the constant Co which normalizes W{x,t) is given by 

^° - 2ViT[a,ay ^^^^ 

and 

p^-^e-Pdp (12) 

is the incomplete Gamma function. 

In the limit that e vanishes, W{x,t) becomes a Gaussian. This can be seen from 

i„n.)~(jij.-i)Mi + .H)-i^~g^H-0(.) (13) 



and hence 



limF(M) ~ exp[-u^/2Do\ (14) 



This is because, in that limit, the diffusion coefficient (191) is a constant. As e increases the 

n 

tails of the distribution decay slower. We refer the reader to Ref. [19|| for a more detailed 
study of the special case e = l/^/Do when F{u) is an exponential distribution. 

We simulated the price returns using random walks with steps of unit size occurring at 



non-constant time intervals. The time between steps is 1/D{x,t) 21[, where D{x,t) was 



calculated at every time step. The simulations consisted of many independent walkers, each 



of which started at the origin, and randomly chose the direction of each event to be either 
an increase or a decrease with equal probability. The walks continued until a maximum time 
was reached. Fig. [T] compares the analytical and simulation results, showing good agreement 
between them. 
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FIG. 1: Distribution of returns for a piecewise linear diffusion coefficient, Eq. Q, with Dq = 1. 
The shape of the distribution changes with the parameter e. Notice the special cases of e = where 
the distribution is Gaussian, and of e = I where the distribution is exponentiab The solid lines 
represent the analytical solution of Eq. (fTU)) . and the data points are the results of the random 
walk simulations with 10^ independent walks, each one lasting a time t = 256. 




B. Quadratic Diffusion 

The second form of D{u) we consider is, 

D{u) = Do{l + eu^). (15) 
It this case, the solution to Eq. ((T)) obtained using Eq. ((Zj), is 

where f3 = l/{2Doe) and the normalization constant for W{x,t) is 

r[i + /?] 



/V^r[i]F[i + /3] 
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This result is plotted in Fig. |2l where it is compared to the results of the corresponding 
discrete random walk simulation for different values of e, with Dq = 1. The simulations were 
performed as in the case of piecewise linear diffusion, except that in this case D is given 
by Eq. (jl5|) . Note also that the results from the simulation are again consistent with the 
analytical solutions. 

As before, the return distribution also becomes a Gaussian in this case when e vanishes. 
However, as e increases, the tails of W{x,t) become power-law distributed. This behavior 
can better be appreciated in Fig. El where a log- log plot for different values of e is presented. 
In the limit of e — > oo the tails of the distribution are well fitted by a power-law with 
exponent 2. Meanwhile, as e ^ the tail can also be fitted with a power-law, but with an 
exponent whose value increases and e decreases. However, the fit is good over a range that 
shrinks as e decreases. This is expected since the distribution becomes Gaussian in the limit 
e — » 0. It is important to point out that as e is decreased from oo to the exponent observed 
in the tail varies from 2 to oo. Thus, these results reproduce the empirical observations of 
real markets that find fat tailed price return distributions with exponents ranging from 2 
upward. Exponents as large as 7.5 have been reported, |^, but at these values the results 
have large error bars. This is because large exponents are found when the time scale is 
increased, and the amount of data samples used in the analysis decreases. 

Fat tails f{x,t) ~ |a;|~"~"^ for |x| >> 1, in the range 2 < a < cxd are also generated by 
symmetric Levy distributions 



exponent H < ^. Here, there is no Langevin equation with in x {D{x,t) doesn't exist), and 
the efficient market hypothesis is violated j^]. 

IV. DYNAMIC CONSENSUS PRICE S 

As has already been mentioned, the price of the stock will fluctuate about the consensus 
value S. However, contrary to what has thus far been considered, there is no reason to expect 
that the consensus value should stay fixed. In fact, it is plausible that with every trade S{t) 
will shift by a small amount toward the value of the current price S{t), or equivalently, that 




(18) 



when (Ax^) = At^/" for a > 2, but this implies fractional Brownian motion with Hurst 
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FIG. 2: Distribution of returns using the quadratic diffusion coefficient in Eq. 1)151) . Note that 
the distribution is not exponential, but instead has fat tails. The solid lines represent the exact 
solution to Eq. (|16|) . and the data points represent the results from the simulation. For e = 0.05 
2 X 10^ walks were simulated, 5 x 10^ for e = 0.5, and 6 x 10^ for e = 1. The final time used in 
each case was t = 256. 
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FIG. 3: Log-log plot of the distribution of x using a quadratic D{x, t) showing that the power- 
law tails of iy(x,t) can have exponents ranging from 2 upward. The data points represent the 
analytical solution, from Eq. 1161 The tail of each case is fitted with a straight line having a slope 
equal to the tail exponent. 

the consensus value of returns x{t) will shift toward the value of the current return x{t). Of 
course, the diffusion constant will change with the consensus value changes to -D(x — x,t) 
in order to keep its minimum at a; = a;. 
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In this section, results of random walk simulations which utilize the non-constant diffusion 
coefficients considered in the previous section, and which allow for the consensus value to 
change, are presented. A very simple dynamics for the value the consensus return x is 
considered; the change in the value at each time step. Ax, is assumed to be proportional to 
the difference in the current return value and its current value. 



where > is a constant that we will assume to be small. There are two essential differences 
between the simulations discussed in section III and those in this section. First, as mentioned 
above the diffusion constant used is D{x — x,t) instead of D{x,t). Therefore, the time 
between steps becomes 1/D{x — x,t). Second, the value of x is varied dynamically using 



The simulations again begin with the consensus price of the stock equal to its initial 
value 5^(0) = Sq, and therefore the initial value of the consensus return vanishes x(0) = 0. 
Subsequently, x will fluctuate around its initial value. Of course, x{t) will also fluctuate 
about the origin. When the value of x{t) is near the origin, it is often the case that |a;(t)| ~ 
|x(t)|. This causes the peak in the distribution of returns W{x,t) to smear out. 

A. Linear Diffusion 

To understand the effects of a dynamic consensus value x on the distribution of price 
returns, first consider the case of piecewise linear diffusion D{x — x,t). As discussed earlier, 
if e = ^= and x is static, this form of the diffusion coefficient will result in an exponential 



return distribution. Figs. |3] and IHl present the results of simulations with dynamic x. As 
expected from the argument in the previous paragraph, the effect of the dynamics of x is to 
smooth out the peak in W{x,t). In fact, it becomes Gaussian in the center, as can be seen 
from the fit to the quadratic function shown in red in Fig. HJ That function is fit through 
the 31 points at the peak of W{x, t). The range of the quadratic region is directly related to 
the value of k. If k increases this region is extended to a larger range, see Fig.O Away from 
the center of the distribution, where the effects of the dynamics of x become less important, 
the exponential form of W{x, t) is retained as expected. This is shown by the fit to the line 
shown in blue, which works in the tail of the distribution. 



Ax(t) = k[x{t) -x{t)], 



(19) 



Eq.dlHD. 
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FIG. 4: Distribution of returns with linear diffusion and dynamic consensus value. The center of 
the distribution of x{t) is well described by a quadratic function. The curve shown was fit to the 
31 points in the center. The tails of the distribution have an exponential decay. A straight line 
was used to fit the tails. The figure shows results for e = 1, Dq = 1 and k = 0.01. 1.6 x 10^ walks, 
each with a final time of t = 256, were used in the simulation. 
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FIG. 5: Dependence on k of the size of the Gaussian region at the center of the distribution. As 
k increases the width of the Gaussian region is increased. In the limit of A; = 1 the distribution 
becomes completely Gaussian. The data points represent the simulation and the solid lines repre- 
sent the quadratic fits. Each simulation consisted of 10'' random walks, with e = 1, Dq = 1 and 
t = 256. For clarity the distributions were shifted vertically and the tails are excluded from the 
plot. 
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Fig. ini shows the distribution of x for different values of /c, which we will call P{x,t) to 
distinguish it from W{x,t). As k is decreased, xq stays closer to the origin, its starting posi- 
tion, and the tails of the distribution decay rapidly. This is why the tails of the distribution 
of X are not affected by the movement x. In the limit of k = 0, x becomes static, and the 
distribution will be a single point at the origin. On the other hand, in the limit k goes to 1 
the distribution becomes a Gaussian. 
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FIG. 6: Distribution of x using the piecewise linear diffusion coefficient. Note that the distribution 
depends on the parameter k. For small k, x becomes more localized around the origin. In this 
simulations, e = 1, Dq = 1 and t = 256 with 10^ walks each. The solid lines are guides to the eye. 

B. Quadratic Diffusion 

Now consider the effects of the dynamics of x on the return distribution for the case where 
D{x — x,t) is a quadratic function. In this case, as we have seen, if x is static, then the 
center of the re turn distribution has a peak, but does not have discontinuity in the slope 
at X = 0. Fig. [7| shows the return distribution calculated from simulations with dynamic 
X. As expected, the peak at x = is broader than in the case of static x, and it can also 
be fitted with a quadratic function, indicating a Gaussian peak. Notice, though, the tail 
behavior in this case differs from that observed when x was static. In this case, the tails of 
the distribution are exponential. We will return to this point at the end of this section. 

Fig. ISl shows the distribution of x. As the value of k decreases, the distribution of x 
becomes sharply peaked. This occurs at the same time that the tails of the distribution are 
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FIG. 7: Distribution of x using a dynamic consensus value and the quadratic diffusion coefficient. 
The central region of can be fitted using a quadratic function (solid line). To obtain the fitted line, 
31 points at the center were used to fit the distribution. The data points show the results from the 
simulation using e = 1, Dq = 1, t = 256, k = 0.0005 and 2 x 10'' walks. 



getting heavier, implying that compared to Fig. IHl there is a larger chance of x being far 
from its original position. This behavior is presumably due to the effect of the fat tails in 
the distribution of x{t) for static x. When the value of x{t) is in the tail of its distribution, x 
is "pulled" far away from the origin. This dynamics is very different than what we observed 
in the case of linear diffusion, where the tails of the distribution of x decayed faster as k 
decreased. Notice that in the limit of A; = the distribution of x will also be a single point 
at the origin, as is also the case for linear diffusion. 

We now take a closer look at the tails of W{x,t). Fig. IHl presents a log-log plot with 
results from simulations for quadratic diffusion using different values of k and e = 1. It is 
observed that as k decreases, the power-law behavior starts to emerge in the tails. In the 
limit of = the results should be the same as in the static x case (a power-law with slope 
3). To explain why the power-law disappear with an increase of the parameter k we turn 
to the dynamics of x{t). When k is increased the value of x will follow closer x{t) making 
D{x — X, t) have a more constant value. This results in the tails of x{t) becoming Gaussian. 
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FIG. 8: Distribution of x using the quadratic diffusion coefficient with dynamic consensus value. 
The tails of the distribution get heavier as k is decreased. This results were obtained using e = 1, 
Dq = 1 and t = 256. The number of walks used was 2 x 10'^ for A; = 5 X 10"^ and k = 10'^, 4 x 10^ 
for A; = 5 X 10^^ and 5 x 10^ for k = 10^^. Here again the solid lines serve as guides to the eye. 

V. CONCLUSIONS 

We have presented a theory for the distribution of stock returns. It is based on the 
conjecture that the rate of trading of a stock depends on how far its current price is from 
a consensus price, S. The resulting models use a non-constant diffusion coefficient D{x,t) 
to simulate the rate of returns. When S is fixed and a piecewise linear coefficient is used, 
an exponential distribution of returns is found. With quadratic diffusion distributions with 
fat tails are found. The exponents describing the power law fat tail distributions range 
from 2 to oo. In both cases we obtained an exact solution for W{x,t) and simulations that 
support our findings. When S is allowed to move, both forms of diffusion coefficient give 
distributions with an approximately Gaussian near the origin. Finally, we note that the 
range of behaviors observed here with this simple model covers the range of non-Gaussian 
behaviors seen in the distribution of returns of real financial markets. 

A.L.A. and G.H.G are partially supported by the Institute for Space Science Operations 
(ISSO) at the University of Houston. K.E.B. is supported by the NSF through grants DMR- 
0406323 and DMR-0427538, and by the Alfred P. Sloan Foundation. G.H.G. is supported 
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FIG. 9: Log-log plot of the distribution of x with quadratic diffusion coefficient and dynamic 
consensus value. The fat tails in the distribution go to a power law as k goes to zero. Each 
simulation used e = 1, Dq = 1, t = 256 and 2 x 10^ random walks. The black line at the top has a 
slope = 3. The distributions have been shifted vertically for clarity. 

by the NSF through grant PHY-0202001. 
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